High precision quantum control of single donor spins in silicon 
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The Stark shift of the hyperfine coupling constant is investigated for a P donor in Si far below 
the ionization regime in the presence of interfaces using Tight-binding and Band Minima Basis 
approaches and compared to the recent precision measurements. The TB electronic structure cal- 
culations included over 3 million atoms. In contrast to previous effective mass based results, the 
quadratic Stark coefficient obtained from both theories agrees closely with the experiments. This 
work represents the most sensitive and precise comparison between theory and experiment for single 
donor spin control. It is also shown that there is a significant linear Stark effect for an impurity near 
the interface, whereas, far from the interface, the quadratic Stark effect dominates. Such precise 
control of single donor spin states is required particularly in quantum computing applications of 
single donor electronics, which forms the driving motivation of this work. 

PACS numbers: 71.55.Cn, 03.67.Lx, 85.35.Gv, 71.70.Ej 



Solid state quantum computer architectures have re- 
ceived considerable attention in recent years due to their 
promise of scalability and utilization of the vast knowl- 
edge and experience of the semiconductor fabrication in- 
dustry. Several proposals employ the electronic states 
of phosphorus donors in silicon to encode qubits based 
on nuclear spin [l[ , electronic spin @, 0, 0] , or the elec- 
tronic states of a singly ionized two-donor system 
Two-dimensional architectures for such qubits have been 
proposed to demonstrate scalability strategies @. 

In Kane's original nuclear spin qubit proposal [H, and 
the more recent semi-global control electron spin pro- 
posal 0], single-qubit operations are performed by elec- 
trically controlling the hyperfine interaction between the 
donor electron spin and the nuclear spin of the P impu- 
rity. To implement rigorous large-scale quantum algo- 
rithms protected by Quantum Error Correction (QEC), 
the donor wave function dependent hyperfine coupling 
constant needs to be controlled to a very high degree of 
precision. Recent experiments have made tremendous in- 
roads towards achieving such control of donor spin. In 
Ref 0] coherent oscillations of P donor electron spins 
were measured by electrically detected magnetic reso- 
nance. Similar techniques were used in Ref [8] to invesit- 
gate the hyperfine interaction in the presence of lattice 
strain for Si grown on Sii-^Ge^. This work is motivated 
by a recent ESR experiment Q that measured the Stark 
shifts for 121 Sb donors in silicon buried at a depth of 
150 nm from an interface. A quadratic dependence of 
the hyperfine coupling constant and electron g-factor on 
the electric field was reported far below the ionization 
regime. A weak linear Stark effect was also reported. In 
comparison with theory, effective mass theory (EMT) in- 
ferred from the work of Friesen [l6[ predicts a quadratic 



Stark shift some ten times greater than the measured re- 
sult. This letter investigates the Stark shift in the experi- 
mentally probed regime well below the ionization thresh- 
old and shows that two sophisticated methods, namely 
Tight-Binding (TB) and Band Minima Basis (BMB), are 
able to predict the system control parameters in contrast 
with EMT. The Sb and P donors in Si can be treated 
as analogous because their ground state binding energies 
differ only by 2.8 meV, and their Is manifold splittings 
due to valley-orbit coupling and local perturbations are 
also similar in nature jl0[ . This is important as both the 
TB and BMB methods employed here are optimized for P 
donors. This agreement between theory and experiment 
for a non-trivial measure of wave function control, by 
two distinct theoretical approaches, represents the most 
precise test of our understanding of donor systems. The 
current level of theoretical understanding of the system 
reported here, in conjunction with the advancement in 
single atom fabrication technologies and the above- 
mentioned experiment, represents one of the first steps 
towards experimental realization of high precision con- 
trol of donor spins in Si. 

Most of the approaches based on EMT have involved 
using hydrogenic envelope functions [12[ with Bohr radii 
calculated analytically or fitted variationally considering 
the anisotropic effective masses of Si. While Kettle et al 
[IK used EMT to study the effects of electric fields and in- 
terfaces non-perturbatively with a TCAD gate potential, 
this work ignored valley-orbit (VO) interaction responsi- 
ble for lifting the degeneracies of the Is manifold. A num- 
ber of authors [3, EH, [l6[ used EMT to investigate the 
donor energy spectrum using symmetry arguments, per- 
turbation theory, single or the multi- valley Schroedinger 
equation. In Ref 0, Friesen considered valley-orbit cou- 
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pling and field effects in the effective mass formalism, to 
solve a multi- valley Schrodinger equation for the six Is 
states in the valley basis. Martins et al [l7| used a TB 
sp 3 s* second nearest neighbor model to study effects of 
electric fields. In Ref [18] the BMB approach was pro- 
posed and applied to the same problem. Calderon et 
al [lH used EMT to investigate shuttling times of the 
donor electron between a donor-bound state and a sur- 
face bound state. Both TB and BMB approaches con- 
sider core effects and VO coupling explicitly to model 
impurities. The TB and the BMB methods are very dif- 
ferent from each other in the sense that TB represents the 
Hamiltonian in real space while BMB represents it in mo- 
mentum space. The close agreement of the results from 
two such different methods of calculations, as shown in 
this work, emphasizes the fundamental nature of the im- 
purity physics. The work presented here resolves signifi- 
cant limitations of the previous treatments by including 
realistically large device components with over 3 million 
atoms in an atomistic treatment under realistic boundary 
conditions. This is also one of the first theoretical works 
on the Stark effect for P impurities to be benchmarked 
with a relevant experiment. While most theoretical ap- 
proaches have concentrated on donor spectrum and wave 
functions near the ionization regime, little work has been 
done to explore the linear and quadratic Stark shifts of 
the hyperfine coupling far below the ionization regime for 
realistically large systems in the presence of interfaces. 

The spin Hamiltonian of a donor electron in an applied 
electric field e and a magnetic field Bqz is, 

H z = g&fi B B S z + A(e)I z S z (1) 

where A{e) and g(e) are the hyperfine coupling constant 
and electron g-factor, respectively, as a function of elec- 
tric field. The first term in the Hamiltonian is the spin- 
orbit Stark shift, while the second term is the hyperfine 
Stark shift. An electric field distorts the shape of the 
donor wave function, and modifies its angular momen- 
tum, which in turn modifies the g-factor. The field also 
pulls the donor wave function away from the impurity 
site, reducing A(e) which is proportional to \^f(e, r~o)| 2 , 
where r~o is the impurity site. The change in A{e) is pa- 
rameterized as, 

AA(s)=A(0)(r ]2 s 2 + r ]l s) (2) 

The quadratic Stark coefficient 772 was measured to be 
—3.7 x 10 -3 /im 2 /V 2 from the experiment by Bradbury 
et al [9]. The EMT result inferred from [16] in Ref @ was 
-2 x 10 -2 /im 2 /V 2 . Since the spin-orbit Stark effect is 
small for practical magnetic fields, this letter focuses on 
the behavior of A(e) only. An impurity placed far away 
from interfaces in an unstrained Si lattice experiences a 
dominant quadratic Stark effect for both A{e) and g(e). 
However, presence of nearby interfaces or lattice strain 
can produce a non- negligible linear Stark effect [13, Hl[ . 



The origin of the linear and quadratic Stark effects can 
be explained qualitatively from perturbation theory. The 
unperturbed ground state wave function for an impurity 
placed far away from interfaces has an approximate even 
symmetry. First order perturbation theory does not yield 
any energy correction if the perturbing Hamiltonian of 
a constant electric field is of odd symmetry. However, 
second order perturbation theory produces a quadratic 
dependence of the corrected energy on the field. If the 
impurity is very close to the interface, the even symmetry 
of the unperturbed wave function is broken. In such a 
case, first order perturbation theory also yields a linear 
dependence of the corrected energy on the electric field. 

TB and BMB methods are separately used to calcu- 
late A{e). From each calculation, the ground state wave 
function ^(s,r) is calculated in the presence of the field 
e. The hyperfine coupling A(e) is determined by, 

m = i*(g>o)i 2 m 

A(0) |¥(0,ro)|2 w 

In the range of low electric fields considered here, the ex- 
cited states (E and T) are not expected to affect A(e) 
since this manifold is separated by 10 meV from the 
ground Ai state. While EMT based approaches are 
concerned with contributions from valley-minima states 
only, TB and BMB approaches consider a more exten- 
sive Bloch structure of the material. TB and BMB also 
include valley-orbit (VO) interaction consistently, while 
most EMT treatments of VO are in violation of some of 
the assumptions made to derive the effective mass equa- 
tion 0. 

The BMB technique solves the Hamiltonian including 
an external electric field in a large but truncated basis of 
pure crystal Bloch states near the conduction band min- 
ima of a host obtained by the pseudo-potential method. 
Since the discretized Schrodinger equation is solved in 
reciprocal space, the Fourier transform of the impurity 
potential is used. Near the impurity core, the potential 
is modified from its Coulombic nature to include central 
cell effects. The corrected potential produces a broader 
Fourier spectrum in k-space and couples different valleys 
to lift the Is degeneracies. The core potential in momen- 
tum (q) space used in [18] is, 

TT ( \ - 1 ( *Cg 2 , k(1 ~ PQ 2 Q 2 n /.x 

Ucor{Q) — o V 9 , 9 ~r 9 , 09 9 1 9) \ ) 

where n is the dielectric constant of Si, assuming a value 
of 11. 9^o- The form and the parameters (£, a, /?, 7) of 
this potential have been taken from Pantelides [22j. An 
overall quenching factor provides a one-parameter fit to 
the ground state energy from which the excited state en- 
ergies and degeneracies compare well with experiment. 
After inclusion of the external field the hyperfine cou- 
pling strength is obtained. 

The standard technique in semi-empirical TB [23[ is 
to optimize a parameter set to accurately reproduce the 
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bulk band structure of a host. Once this is done, the 
same set is used for atomistic modeling of any devices 
made of that host. For this work, the TB parameter set of 
sp 3 d 5 s* nearest-neighbor spin model for Si was optimized 
by genetic algorithm [24]. The donor potential used is, 



V(r) 

47Tft|r — r~o| 

V(r) = U , 



r = f 



(5) 
(6) 



The on-site core correcting potential Uq is adjusted and 
found to be 4.33 eV to obtain the experimental ground 
state energy of 45.6 meV for Si:P. The potential due to 
the electric field is added to the diagonal of the Hamilto- 
nian. The real space Hamiltonian with closed boundary 
conditions is solved by a parallel Lanczos algorithm in 
NEMO 3D for the donor wave function. 

Both methods described here are applicable to other 
hosts and impurities. To describe other Group V impu- 
rities, the hydrogenic coulomb potential is still valid in 
the bulk, but the core corrections need to be adjusted to 
reflect correct experimental binding energies. 
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FIG. 1: (a) Electric field response of hyperfine coupling at 
various impurity depths (BMB and TB). (b) Quadratic (LH 
axis) and linear (RH axis) Stark coefficients with depth (TB). 
(c) Mean position of the ground state electron distribution 
(dipole moment) as a function of the electric field (TB). (d) 
The electric field gradient of the dipole moments (i.e. the 
slopes of the lines in (c)) with respect to depth (TB). 

Fig. 1 summarizes the effects of the electric field and 
the interface on the donor electron. The TB calculations 
use a domain of 32 nm x 65 nm x 32 nm Zincblende 
lattice with 3.45 million atoms. The distance between 
the impurity and the interface is varied parallel to the 
electric field. The BMB calculation assumes a maximal 
depth of 10.86 nm from the interface, and employs a basis 
set of 7986 states. The choice of this depth is dictated by 
the ease of computation as it is sufficiently deep to nullify 



surface effects while not too deep to make the problem 
computationally intractable. In TB, a range of depths 
from 5 nm to 32 nm have been considered. For each 
TB data point, the typical computation times require 
about 7 hours on 20 CPUs 0. Fig. 1(a) shows the 
variation of /S.A(s) with electric field for various impurity 
depths. The data are fitted to the quadratic equation of 
([2]) . As the depth increases, the quadratic coefficient 772 
approaches a constant value, while the linear coefficient 
771 becomes negligible (Fig. 1(b)). For small impurity 
depths, 771 is comparable to 772, which results in a shift of 
the peak of the parabola in Fig. 1(a) towards a non-zero 
electric field. If the linear Stark effect is negligible, an 
applied electric field has two effects on the ground state 
wave function: 1) a decrease in the peak amplitude of the 
wave function at the impurity site, reflected by a decrease 
in A(e) in Fig. 1(a) for higher depths, 2) a shift in the 
mean position of the wave function opposite the electric 
field, giving rise to a non-zero dipole moment as shown 
in Fig. 1(c). The dipole moments vary linearly with the 
electric field, and their slopes approach a constant value 
as the depth increases (Fig. 1(d)). The wave function 
plots of Fig. 2 also demonstrate these effects. The larger 
distortions of the wave function at higher electric fields 
explain the increasing dipole moments. 
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FIG. 2: (a) Electric field induced differential map of donor 
electron wave probability density, |^(e,^o)| 2 — |^(0,^o)| 2 , 
shown as a 2D cut through the impurity center at zo — 16.29 
nm for e — 0.5 MV/m. The electric field is directed along the 
negative y axis, (b) ID cut though the center of the impurity 
parallel to the electric field showing the differential map of 
the probability density for two different electric fields. 



If the impurity is close to the interface, the wave func- 
tion is distorted even at zero electric field. This distortion 
comes from changes in the slope of the part of the wave 
function between the impurity and the interface. Also 
the mean position of the wavefunction shifts farther as 
the distance between the impurity and the interface is 
reduced. These effects give rise to an asymmetric charge 
distribution about the impurity, and a non-zero dipole 
moment is observed even at e = 0. This is verified by the 
shift of the y-intercept of the lines in Fig. 1(c) for small 
impurity depths. In effect, the interface behaves like an 
electric field pushing the donor electron away from it. 
The strength of this field increases as the impurity is 
placed closer to the interface. This means that a higher 



electric field directed away from the interface is needed 
to counteract these interface effects and to restore the de- 
creasing behavior of A(e), as demonstrated in Fig. 1(a). 

Some of the numerical results can be explained quali- 
tatively with the perturbation model using a basis of un- 
perturbed impurity wave functions {\I/o, ■■■}> 
where the superscript represents the order of correction 
and subscript the wave function index in order of de- 
creasing binding energy (0 being the ground state). For 
an electric field Hamiltoninan given by H = qsy, the 
wave function corrected to first order is expressed as, 



TABLE I: Comparison of the quadratic Stark coefficients from 
experiment, EMT, BMB and TB. 
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where i is the degeneracy index for a state m. The dipole 
moment D = q{^l\y\^Q) is then given by, 
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The quadratic term in e can be neglected compared to the 
other two terms for small e. This results in a linear func- 
tion with intercept and slope depending on (^qM^o) anc ^ 
K*oM*m,;)| 2 (£o - ^m,;)" 1 respectively. For an impu- 
rity far away from the interface, the unperturbed ground 
state wave function is of even symmetry, and the expres- 
sion for the intercept evaluates to zero. Since the un- 
perturbed wave function symmetries and eigen-energies 
are not affected by depth any more, the expression for 
the slope assumes a constant value. On the other hand, 
for small impurity depths, the unperturbed wave func- 
tions get distorted depending on the depth. Since the 
ground state has no longer an even symmetry, the in- 
tercept evaluates to a non-zero quantity, in consistence 
with Fig. 1(c). The transition probabilities between the 
ground state and the excited states are also changed due 
to these distortions, causing the slope of the dipole mo- 
ment to vary with depth, in consistence with Fig. 1(d). 
The same reasoning explains why the quadratic Stark co- 
efficient approaches a constant value and the linear Stark 
coefficient becomes negligible with increasing depth. 

Table 1 gives the comparison between theory and ex- 
periment and shows that the EMT estimate of 772 differs 
by an order of magnitude, while TB and BMB estimates 
agree closely with each other and the experiment. The 
BMB estimate for a depth of 10.86 nm is already close to 
the converged TB estimate, although the peak AA(e) at 
a non-zero e in Fig. 1(a) indicates that the BMB estimate 
would improve even more at greater depths. 

In summary, the Stark shift of the hyperfine coupling 
for Si:P in the presence of interfaces is analyzed from 
two very different theories (BMB and TB). The results 
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150 
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EMT (P) [16] 


00 
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-2.74 x 10" 3 
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10.86 


-2.57 x 10" 3 




21.72 


-2.76 x 10" 3 



of both methods are consistent and agree well with the 
experimental data for Si:Sb, thereby providing the most 
sensitive test of our understanding of shallow donor quan- 
tum control to date. For future work, both TB and BMB 
can be optimized for Sb donors, and the Stark shifts can 
be investigated in the presence of interfaces as well as lat- 
tice strain and a gate structure more closely resembling 
that of the experiment. In conclusion, we emphasize the 
need for very close interaction between theory and ex- 
periment to accomplish high precision control required 
in quantum computing. 
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